PACIAE Model Predictions for Pb + Pb Collisions at LHC Compared to the Au + Au 

Collisions at RHIC 



Ben-Hao Sa 1,2,3 *, Dai-Mei Zhou 2 , Bao-Guo Dong 1 , Yu-Liang Yan 1 , Hai-Liang Ma 1 , and Xiao-Mei Li 1 
1 China Institute of Atomic Energy, P. O. Box 275 (18), Beijing, 1024-13 China 
2 Institute of Particle Physics, Huazhong Normal University, Wuhan, 430079 China 
3 CCAST (World Lab.), P. O. Box 8730 Beijing, 100080 China 

The parton and hadron cascade model PACIAE is used to simulate the 0-6, 15-25, and 35-45% 
most central Au + Au collisions at y / sjvjv = 19.6, 62.4, 130, and 200 GeV and the 0-10% most central 
Pb + Pb collisions at y / 5jvjv=5500 GeV. The calculated charged multiplicity and the charged particle 
transverse momentum distribution, pseudorapidity distribution, and extended longitudinal scaling 
for Au + Au collisions well agree with the corresponding PHOBOS data. Thus the above observables 
calculated for Pb + Pb collisions would be a reliable prediction. 
PACS numbers: 25.75.Dw, 24.85. +p 



In Ref. [l| the predictions of variety models for Pb+Pb 
collisions at LHC energy (y / 7/vw=5500 GeV) have been 
compiled including the thermal (statistical) models 
The basic assumption in thermal model is that the fi- 
nal state hadrons in relativistic heavy-ion collisions are 
originated from a thermal source of a given thermal pa- 
rameters of temperature and baryonic chemical potential 
Q . In 0] a parametrization of the temperature and bary- 
onic chemical potential as a function of ^/snn (in GeV) 
is obtained from the fits of calculated particle yields and 
ratios to the experimental data at mid-rapidity over a 
broad energy range of y / ijvjv= 2. 7-200 GeV. Then they 
extend above parametrization to the LHC energy and 
provide a quantitative prediction for LHC Pb + Pb ex- 
periments. 

Similarly, the transport models are also unable to be in 
agreement with the experimental data of nucleus- nucleus 
collisions in a broad energy range from AGS to SPS, 
RHIC, and even to LHC without adjusting any parame- 
ters. Of course, the less number of adjusting parameters 
the better the model is. In this letter the parton and 
hadron cascade model, PACIAE, is used to simulate the 
0-6% (15-25 and 35-45%) most central Au + Au colli- 
sions at y/sj^= 19.6, 62.4, 130, and 200 GeV and the 0- 
10% most central Pb+Pb collisions at ^/sjvat=5500 GeV. 
The charged multiplicity and the charged particle trans- 
verse momentum distribution, pseudorapidity distribu- 
tion, and extended longitudinal scaling are calculated. 
In the calculations we fix all model parameters, except 
the parameter b in Lund string fragmentation function 

1 

f(z) cx z _1 (l - z) a exp{-bm 2 T /z). (1) 

This function expresses the probability that a given z is 
picked, here z refers to the fraction of string energy (mo- 
mentum) taken away by the produced hadron. In Eq. 
([T|) the rriT = \/ v\ + m? is transverse mass of hadron, 
Pt and are, respectively, the transverse momentum 



and rest mass of hadron, and the a and b are parame- 
ters. As mentioned in [(J that the parameter b is varied 
with string density and the higher reaction energy (tem- 
perature) is corresponding to the lager b. According to 
the experimental facts that the temperature of fireball, 
in relativistic heavy-ion collisions, as a function of ^/snn 
increases dramatically first and then approaching satu- 
ration gradually p}, we assume b to be approximately a 
function of a/sjvtv as 



fe = 6 °^7r< l f V^< 200 GeV, 



b = bo[l + (l --==)], *f V^v"> 200 GeV, (2) 

where &o=6- As the calculated above observables in Au + 
Au collisions are in good agreement with PHOBOS data 
(see later), the calculated those observables in Pb + Pb 
collisions are thus a reliable prediction. 



TABLE I: Charged multiplicity in Au+Au collisions at ^/sjvjv 
=19.6, 62.4, 130, and 200 GeV and in Pb + Pb collisions at 
V5jvjv=5500 GeV. 

Au+Au Pb+Pb 

Energy 19.6 62.4 130 200 5500 

(GeV) 

Ccntrality 0-6% 0-10% 

Multiplicity 

Exp. 1689+100 1 2845 + 142 2 4170+210 1 5060+250 1 
PACIAE 1533 2919 4140 5001 14695 
^+\ v= o 250 3 425 3 559 3 641 3 1196 4 
b 058 2 4 6 12 

1 PHOBOS data taken from [11]. 

2 PHOBOS data taken from [121 . 

3 Estimated from calculations in 0.2< rj <1.4 . 

4 Estimated from calculations in |?7| <0.12 . 
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The parton and hadron cascade model, PACIAE [8|, is 
based on PYTHIA (the model for hadron-hadron colli- 



2 



10 
10 

cf lO 1 

I 10 

o io 
£ io 

•a 

f io- 
S l0 + 

J 3 -4 ! 

U 10 r 
n 10 r 

S -f 

10 
10"' 

io- ! 



Au+Au @ /s NN =200GeV^ 

(h + +h)/2 
0.2<T1<1.4 



■ i< :v.,'iii i 
if 



(a) z j 

v ; A35-45%(»10 ) 

Full labels: PACIAE A 
Open labels: PHOBOS A 



Au+Au 8 /s^ N =200 GeV 
Au+Au @ /s^,=200 GeV 
Pb+Pb @ /sT N =5.5 TeV 
(h + +h)/2 



T 0-10« 

▼ ▼ 



(b) 

■ 

full Tj phase space 



reads 



l 



2 3 4 
p [GeV/cl 



5 



6 8 10 12 14 16 
P T [GeV/cl 



FIG. 1: The charged particle transverse momentum distribu- 
tions in Au + Au collisions at y / 5jvjv=200 GeV and Pb + Pb 
collisions at x /sjvjv=5500 GeV. PHOBOS data are taken from 
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where the factor K is introduced counting for the 
higher order pQCD and non-perturbative QCD cor- 
rections, a s stands for the strong (running) cou- 
pling constant, and s, t, and u are the Mandelstam 
variables. For the process q\qi — » 9192, for instance, 
one has 
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(4) 



91 927+9192 

It diverges at t—0 and has to be regularized by 
introducing the parton colour screen mass \i as fol- 
lows 



sion) 0] and is composed of four stages: the parton ini- 
tialization, parton evolution (scattering), hadronization, 
and hadron evolution (rescattering). 

1. Parton initialization: 

In the PACIAE model the nucleus-nucleus collision 
is decomposed into nucleoli -nucleon (NN) colli- 
sions according to the geometry of nucleus- nucleus 
collision. A NN collision is performed by the 
PYTHIA model [5| with the hadronization process 
switched-off. Thus the consequence of a nucleus- 
nucleus collision is a state composed of quark 
pairs, diquark pairs, gluons, and very few hadronic 
remnants. If the diquark (anti-diquark) is split 
forcedly into quarks (anti-quarks) randomly, the 
partonic initial state of a nucleus- nucleus collision 
is reached. 



2. Parton evolution (scattering): 

The parton initialization stage is followed by parton 
evolution (scattering). Here the 2 — » 2 LO-pQCD 
differential cross sections @ are employed. The 
differential cross section of a sub-process ij — > kl 
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(5) 



The total cross section of the parton collision i + j 
is then 



kj J ~ S 



da. 



ij — >kl 



dt 



(6) 



With above total and differential cross sections the 
parton evolution (parton scattering) can be simu- 
lated by the Monte Carlo method. 

Hadronization: 

The hadronization at the moment of partonic 
freeze-out (no more partonic collision at all) is con- 
sequent on the parton evolution stage. In the PA- 
CIAE model, the phenomenological fragmentation 
model (String Fragmentation, Independent Frag- 
mentation, or Cluster Fragmentation) || and the 
coalescence model are supplied for hadronization of 
partons after scattering. The String Fragmentation 
model is adopted in this letter. We refer to [|| for 
the details of the hadronization stage. 



4. Hadron evolution (rescattering): 

After hadronization the rescattering among pro- 
duced hadrons is dealt with the usual two-body col- 
lision model. The details of hadronic rescattering 
can see [lOj . 

We compare the calculated charged multiplicity with 
the PHOBOS data (taken from [II El) in 0-6% most 
central Au + Au collisions at y / Sjv~iv =19.6, 62.4, 130, and 
200 GeV in Table D One sees here that the PHOBOS 



data are well reproduced within error bars by the PA- 
CIAE calculations, except the case of y / sj v lv=19.6 GeV. 
That is consistent with the fact that the PYTHIA model 
is more suitable for higher reaction energy. The predicted 
charged multiplicity and charged particle pseudorapid- 
ity density at mid-pseudorapidity in 0-10% most central 
Pb+Pb collisions at y / s]vlv=5500 GeV is also given in this 
table. That pseudorapidity density, ^+y-\ v=a —1200, is 
within the predicted values of other sixteen models listed 
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FIG. 2: The charged particle pseudo-rapidity distributions in 0-6% most central Au + Au collisions at y / ijvjv=19.6, 62.4, 130, 
and 200 GeV and in 0-10% most central Pb + Pb collisions at - v /sjvjv=5500 GeV. The full labels are calculated results by the 
PACIAE model and the open labels are the corresponding PHOBOS data. PHOBOS data are taken from [ll[ (for y / sjvjv = 19.6, 
130, and 200 GeV) and pj (for ^/s^ =62.4 GeV). 




In Fig. Q] (a) we compare the calculated charged par- 
ticle transverse momentum distributions with the cor- 
responding PHOBOS data (taken from [l4|) in Au + 
Au collisions at % /sjv r /v=200 GeV. Here the open cir- 
cles, squares, and triangles are, respectively, the PHO- 
BOS data of the 0-6, 15-25, and 35-45% most central 
Au + Au collisions and the full labels are the PACIAE 
model results. We see in panel (a) that the PHOBOS 
data are reasonably good reproduced. The prediction 
for charged particle transverse momentum distribution 
(triangles-down) in 0-10% most central Pb + Pb colli- 
sions at y / 7/\w=5500 GeV are given in Fig. Q] (b). For 
comparison, the charged particle transverse momentum 
distributions in 0-6 (circles) and 15-25% (squares) most 
central Au + Au at GeV are also given in 

panel (b) . Note that the charged particle transverse mo- 
mentum distributions shown in Fig. [T] (b) is integrated 
over full pseudorapidity. 

Figure[2]shows the charged particle pseudorapidity dis- 
tributions in 0-6% most central Au + Au collisions at 
v / saT/v=19.6 and 62.4 GeV (panel (a)) and at y/sNN=130 
and 200 GeV (panel (b)) and in 0-10% most central 
Pb + Pb collisions at y / s]v]v=5500 GeV (panel (c)). In 
this figure the open labels are the PHOBOS data and the 
full labels are the PACIAE model results. The squares 
and triangles-down in panel (a) are, respectively, for 
v /J/vtv=62.4 and 19.6 GeV and the circles and triangles- 
up in panel (b) are, respectively, for y/sNN =200 and 130 
GeV. One sees again in panels (a) and (b) that the PHO- 
BOS data are reasonably good reproduced, except the 
case of ^5^=19.6 GeV. Fig. H (c) gives the PACIAE 
model prediction for the charged particle pseudorapidity 
distribution in 0-10% most central Pb + Pb collisions at 
Y / s]v7V=5500 GeV. In panel (c) one sees that there is a 
deep valley at the mid-pseudorapidity. 

In Fig. [3] we give the shifted charged particle pseudo- 
rapidity distributions in 0-6% most central Au + Au col- 
lisions at y/sNN— 19.6 (triangles-down), 62.4 (squares), 



130 (triangles-up), and 200 GeV (circles) and in 0-10% 
most central Pb + Pb collisions at y/sNN—5500 GeV (di- 
amonds). Here the shifted pseudorapidity is 

V' = V - Vbeam, (7) 

where ybeam is the beam rapidity. We see in Fig. [3] 
that the extended longitudinal scaling [l2j is well kept 
not only among Au + Au collisions at a variety values 
of y/SNN but also among Au + Au and Pb + Pb collisions. 
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FIG. 3: The shifted charged particle pseudorapidity distribu- 
tions in 0-6% most central Au + Au collisions at % /sjvjv=19.6, 
62.4, 130, and 200 GeV and in 0-10% most central Pb + Pb 
collisions at y / ijvjv=5500 GeV. The shifted pseudorapidity 
rj — r) — ybeam, where yteam is the beam rapidity. And the 
Np ar t refers to the number of participant nucleons. 

In summary, we have used the parton and hadron 
cascade model PACIAE to simulate the 0-6, 15-25, and 
35-45% most central Au + Au collisions at y / i]v r /v=19.6, 
62.4, 130, and 200 GeV and the 0-10% most central 
Pb + Pb collisions at ^/sjv at=5500 GeV. The charged 
multiplicity and the charged particle transverse mo- 
mentum distribution, pseudorapidity distribution, and 
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extended longitudinal scaling are calculated. For 
Au + Au collisions the calculated above observables are 
in good agreement with the corresponding PHOBOS 
data. Thus the above observables calculated for Pb + Pb 



collisions would be a reliable prediction. 
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